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The existence of multiple energy scales is regarded as a signature of the Kondo breakdown 
mechanism for explaining the quantum critical behavior of certain heavy fermion compounds, like 
YbRh2Si2. The nature of the intermediate state between the heavy Fermi liquid and the quantum 
critical region, however, remains elusive. In this study we suggest an incoherent heavy-fermion 
scenario, where inelastic scattering with novel soft modes of the dynamical exponent z = 3 gives 
rise to non-Fermi liquid physics for thermodynamics and transport despite the formation of the 
heavy-fermion band. We discuss a crossover from z = 3 to z — 1 for quantum phase fluctuations. 
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Research on quantum criticality has been one driv- 
ing force in modern condensed matter physics, where the 
universal scaling reflects the non-perturbative nature of 
strong correlations [l|, The observation of a regime 
with T-linear resistivity is the hallmark of quantum crit- 
icality in heavy fermions [sj . This observation combined 
with the presence of anomalous exponents calls for an 
interacting nature of the non- Fermi liquid fixed point [1] . 
Heavy-fermion quantum criticality has been regarded as 
a rule model, where competition between the Kondo ef- 
fect and Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action gives rise to a quantum phase transition from an 
antiferromagnetic metal to a heavy-fermion Fcrmi-liquid. 
Non-Fermi liquid physics is displayed in the quantum 
critical region 0, Q . 

Two competing theoretical frameworks have emerged, 
referred to as the Kondo breakdown (KB) mechanisrn L 
[To| and the spin-density-wave (SDW) scenario [ll|, 
respectively. Although these competing scenarios cover 
the T-linear transport [1, [T^, only the KB mecha- 
nism could explain the divergent Griineisen ratio with 
the special critical exponent of 2/3 in YbRh2Si2 In 
addition, another KB scenario based on the slave-fermion 
representation uncovered the diverging uniform spin sus- 
ceptibility with an exponent 2/3, consistent with an ex- 
periment for YbRhaSia [H. 

As shown in the above discussion, critical exponents 
can be thought as a fingerprint of each scenario. These 
exponents can be found from the Eliashberg approxi- 
mation, where self-energy corrections for both electrons 
and critical fluctuations are introduced self-consistently 
[Tsl [l6| . However, the stability of the Eliashberg frame- 
work has been questioned recently because electrons turn 
out to be strongly interacting at quantum critical points 
(QCPs) even in the large- iV limit, the cornerstone of the 
Eliashberg theory, where N is the number of fermion col- 
ors [13, UM- In ttiis respect it is desirable to flnd non- 
perturbative features beyond the Eliashberg approxima- 



tion. 

Recently, two of us predicted violation of the 
Wiedemann-Franz law at the KB QCP, where the ex- 
istence of additional entropy carriers, identified with 
charge-neutral spinon excitations, gives rise to additive 
contributions for the thermal conductivity, resulting in 
enhancement of the Lorentz number [l9| . Furthermore, 
the KB theory was claimed to show an abrupt collapse in 
the Seebeck coefficient from the KB QCP to the SDW or 
spin liquid phase because breakdown of the Kondo effect 
prohibits spinous from carrying electric currents below a 
characteristic energy scale E*, where Fermi-surface fluc- 
tuations start to be frozen and electrons in the f-orbital 
become localized [l^. These two features are based on 
reconstruction of the Fermi surface at the QCP, distin- 
guishing the KB scenario from the SDW theory undebat- 
ably. 

In this letter we investigate another signature of the 
KB mechanism. The Hall coefficient has revealed a novel 
energy scale T* higher than the Fermi-liquid tempera- 
ture TpL, observed in the heavy-fermion side [2ll,[2§|- It 
seems to show an abrupt decrease at T* , but the non- 
Fermi liquid transport and thermodynamics are still ob- 
served in TpL < T < T*. The abrupt change of the 
HaU coefficient is believed to originate from the Fermi- 
surface reconstruction, and has been corroborated by ob- 
servations of a change in the magnetoresistance and field 
dependence of the magnetization [2^ . 

Introducing the phase variable of the hybridization 
order parameter into the KB theory, we propose that 
the intermediate region is characterized by an incoherent 
heavy-fermion band, where quantum phase fluctuations 
give rise to incoherent scattering of heavy electrons and 
do not allow their Fermi liquid behaviors. This preformed 
heavy-fermion scenario shows similarities with the pre- 
formed pair scenario for the pseudogap phase of high Tc 
cuprates [2^ . 

We start to discuss the Kondo effect in the single impu- 
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rity problem. As well known, the slave-boson mean-field 
theory allows for a strong coupling fixed point, identified 
with the local Fermi liquid state [23] ■ However, it causes 
an artificial second order transition at finite tempera- 
tures, which should not exist in the single impurity prob- 
lem. Fluctuation corrections are introduced to check the 
stability of the mean-field state, where they can be iden- 
tified with contributions from vertex corrections to the 
boson condensation [2^. Such soft modes cause the in- 
frared Zog- divergence, argued to make condensation pro- 
hibited, where an infinite-order summation based on the 
parquet approximation will turn the Zog-divergence into 
a power-law behavior. On the other hand, this treatment 
turns out to recover correlation functions such as the spe- 
cific heat coefficient and spin susceptibility of the local 
Fermi liquid. 

We apply this scheme to the heavy-fermion problem, 
described by an effective Anderson lattice model 

= H 4aidr - ^i)c^a - tY^icla^ja + H.C.) 

+Vj2{dlo.a + H.C.) + J2 dlidr + ef)d,, 

i i 

+jJ2s^-s,, (1) 

which shows competition between the Kondo effect (V) 
and the RKKY interaction (J). Cia represents an elec- 
tron in the conduction band with its chemical potential 
fi and hopping integral t. dia denotes an electron in 
the localized orbital with an energy level e/. The lo- 
calized orbital experiences strong repulsive interactions, 
thus either spin-f or spin- J, electrons can be occupied at 
most. This constraint is incorporated in the U(l) slave- 
boson representation, where the localized electron is de- 
composed into holon and spinon, dia = b\fia, supported 

by the single-occupancy constraint b\bi + f-^fia = SN in 
order to preserve the physical space. S" = 1/2 is the size 
of spin and N is the spin degeneracy, where the physical 
case is iV = 2. 

Rewriting the Anderson lattice model in terms of 
holons and spinous, we obtain 

+^ E(^'/-^- + ^•^■) + E ^J^-^i 

^ i i 

+ E fli^r + - JE(/-X.,7,. + H.C.) 

+iJ2x,(blb, + flu - SN) + |x»,f , (2) 

where the RKKY spin-exchange term for the localized or- 
bital is decomposed with the single occupancy constraint 



via exchange hopping processes of spinous with a hop- 
ping parameter Xij , and is a Lagrange multiplier field 
to impose the single-occupancy constraint [27]. 

The saddle-point analysis with bi — >■ b, Xij ~^ X: and 
«Ai — > A reveals a breakdown of the Kondo effect, where 
a spin-liquid Mott insulator (6 = 0) arises with a small 
area of the Fermi surface in J > Tk while a heavy Fermi 
liquid (6 7^ 0) obtains with a large Fermi surface in Tk > 

J [S-Q. Here, Tk — D exp ^ pfp^yi ^ is the single-ion 

Kondo temperature, where pc ~ (2D)~^ is the density of 
states for conduction electrons with the half bandwidth 
D. Reconstruction of the Fermi surface occurs at J ~ 
Tk- 

Quantum critical physics is characterized by critical 
fluctuations of the hybridization order parameter, intro- 
duced in the Eliashberg theory 28] . Dynamics of critical 
Kondo fluctuations is described by z = 3 critical theory 
due to Landau damping of electron-spinon polarization 
above an intrinsic energy scale E* while by z = 2 di- 
lute Bose gas model below E* [1, Q. Here, z is the dy- 
namical critical exponent, which tells the dispersion of 
bosonic modes. The energy scale E* originates from the 
mismatch of Fermi surfaces of conduction electrons and 
spinous, one of the central aspects in the KB scenario. 
Physically, one may understand that quantum fluctu- 
ations of the Fermi-surface reconfiguration start to be 
frozen aX T k, E* , thus the conduction electron's Fermi 
surface dynamically decouples from the spinon's one be- 
low E*. 

We point out that the mean-field transition from the 
z = 3 quantum critical region to the heavy-fermion phase 
is identified with T* of the Hall coefficient [21I, [23| , where 
quantum phase fluctuations of the holon order parameter 
reduce the Fermi liquid temperature much. Decompos- 
ing the hybridization order parameter into its amplitude 
and phase, b = yTp^e^^'' , and performing the continuum 
approximation in terms of low energy fluctuations, we 
reach the following expression 

£ = Cl{dr - Pc)Ca + T: \{dr + iAr)Ca\'^ 

Zrrir 



Me"''c*f, + H.c.) 



V 
7^ 



+f*{dr - + e/ + A - iar)fa + -z \{dr ~ iar)!aY 

2m/ 

+ipb{drOb — lA — ar) + - — [drOb — iX — ar)'^ 
2ub 

Pb 



+- — (drOb-ar-Ar) + —^{di_,a^ - O^a^) 



2mb 



45^ 



(3) 



mc = ^ is the band mass of conduction electrons, and 



Ar is an electromagnetic field, to/ = is the band 
mass of spinous, and A is the mean-field value of the La- 
grange multiplier field with its fluctuation part a-r. ar 
originates from the angular part of the hopping param- 
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eter, Xij = xe^^'^ j playing the role of the U(l) gauge 
field, mb ~ yrjN band mass of holons, origi- 

nating from the electron-spinon polarization function at 
high energies. The low energy physics from such Fermi- 
surface fluctuations is given by the Landau damping term 
in the holon (phase-fluctuation) propagator [Eq. (7)]. 
Uh is a coupling constant for local interactions between 
holons, given by and phenomenologically intro- 

duced, and the second-order time-derivative term with 
Ub results from integration of Spb with pb — > Pb + Spt 
[13] ■ g is the gauge-matter coupling constant. 



This effective field theory is reduced to the slave-boson 
mean-field theory when phase fiuctuations are neglected, 
where pb is identified with 6^. Thus, the mean-field tran- 
sition temperature is identified with T* because a fi- 
nite value of Pb generates the heavy-fermion band, the 
Hall coefficient being reduced due to the Fermi-surface 
reconstruction. If e**'' « I + iOb is performed in the 
Kondo-interaction term of the single impurity problem 
and phase fiuctuations are integrated over up to the sec- 
ond order, we can see that an additional Zog-divergence in 
the spinon self-energy cancels the Zog-divergence in the 
holon condensation, allowing the amplitude (pb) of the 
holon condensation to be finite [131 . This means that the 
mechanism for disappearance of the holon condensation 
lies in transverse (phase) fiuctuations [1^. On the other 
hand, such Goldstone modes turn out to be not harm- 
ful for ordering in the heavy-fermion problem with three 
dimensions. As a result, the heavy-fermion Fermi-liquid 
state is stable against gaussian fluctuations of Goldstone 
bosons 9b. However, the stability is not guaranteed any 
more if quantum phase fiuctuations are taken into ac- 
count beyond the gaussian order. The non-linear a model 
approach is convenient to describe interactions between 
phase modes [1^ , where the phase factor is replaced with 
a complex variable V'- This complex field should be con- 
strained with the unimodular condition, — /i^d-f/jp — 1) 
introduced into the effective Lagrangian, where p,^ is an 
effective chemical potential. 

Rewriting the effective Lagrangian Eq. (3) in terms of 
ip, and introducing quantum corrections self-consistently 
in the Luttinger-Ward functional approach [l^, we ob- 
tain coupled equations for self-energy corrections of elec- 
trons, spinous, phases, and gauge fields. Since vertex cor- 
rections are not taken into account, these self-consistent 
equations are essentially the same as those of the quan- 
tum critical regime in the KB theory [1, A novel 
feature beyond the previous consideration is to introduce 
an additional energy scale p^ , describing coherence of the 
heavy-fermion band. The formation of the heavy-fermion 
band is determined from pb, controlled by A. 

We derive self-consistent equations for three order pa- 
rameters from the Luttinger-Ward free energy functional 
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FIG. 1: A phase diagram in the preformed heavy-fermion 
scenario, where QC, HF, and FL denote quantum critical, 
heavy fermion, and Fermi liquid, respectively. T*{V) corre- 
sponds to the mean-field transition temperature (pf, — 0) in 
the Kondo breakdown theory while the Fermi-liquid temper- 
ature Tfl{V) is much reduced due to quantum phase fluctua- 
tions of the hybridization order parameter ((e**'') = 0). Reen- 
trant behaviors are found in both T* and Tfl numerically, 
but it is not clear whether this effect is fundamental or not due 
to quantum fluctuations. pfff{V,T*) and pbiV^Tph) are also 



shown, where ttIc 
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and blue-circle lines and = Ai^ = 50m ^ for the green- 
square line [27|. The unit of each ajcis is mj ^. 
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Gf{k,iuj) = NS, 



(5) 



where G f(^^j{k,iuj) is the renormalized Green's function 
of spinous (phases) with the heavy-fermion band and 
gc{k,iuj) is the bare Green's function of electrons. V^^''^^ 

and K-'^p'^^ are renormalized Fermi velocity and renor- 
malized Fermi momentum of electrons (spinous) in the 
heavy-fermion band, respectively. p^j:j — p^ — pbX — 
S^(0, 0) is an effective chemical potential, which de- 
termines the Fermi-liquid temperature TpL, where the 
constant contribution of the ip self-energy is I]^(0,0) = 

We perform the numerical analysis, where self-energy 



4 



corrections are evaluated analytically. A detailed proce- 
dure can be found in our supplementary material [27| . 
Figure 1 displays the intermediate state, where pb is 
finite, resulting in the formation of the heavy-fermion 
band, while its coherence is not achieved yet, reflected in 
the fact that the chemical potential —p-'t.ff ^ ^- 
characterized by Pb{T*) — 0, and TpL is determined by 
4/(T^l)=0. 

In the preformed heavy-fermion state the self-energy 
correction of ij) is governed by Landau damping from in- 
coherent heavy fermions, AE^(g,ir2) — 7-^, where the 

damping coefficient is given by 7 = pb ,J,.e ■ Then, 
the imaginary part of the propagator becomes 

-<iG^{q,n) = ^ (7) 

This expression displays a crossover from z — 2> to z = \ 
at Ti sa i \/^^^~f^eff)^ ^1 remains larger than 

Tfl- In Tfl <T < Ti, it is given by -QG^{q,n) = 

Inserting the z — I propagator into self-energy equa- 
tions for fermions, one finds that scattering with such 
fluctuations is less relevant for self-energy corrections of 
fermions than Fermi-liquid corrections in three dimen- 



sions. As a result, we expect that the T-linear resistivity 
due to scattering with z = 3 critical modes ^] becomes 
smoothly transformed into the Fermi-liquid resistivity in 
the intermediate phase. 

Recently, the T* line was proposed to be a Lifshitz 
transition [30|. motivated by the observation that iso- 
electronic chemical doping does not change T* while 
it affects the Neel temperature seriously [2j|- On the 
other hand, non-isoelectronic chemical doping changes 
T* clearly, when Rh is replaced with Fe [3l| . We believe 
that this issue should be clarified. 

In conclusion, we uncovered a new incoherent heavy- 
fermion state which can be relevant to the nature of the 
intermediate region of Tp^ < T < T* . The mecha- 
nism turns out to be existence of quantum phase fluc- 
tuations in the hybridization order parameter. Despite 
the formation of the heavy-fermion band, this interme- 
diate state will show non-Fcrmi liquid physics in trans- 
port and thermodynamics due to scattering with such 
z — i soft modes. The non- Fermi liquid physics become 
transformed into the Fermi liquid physics continuously, 
as the 2 = 3 critical mode turns into z — I, irrelevant for 
fermion dynamics. 
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Appendix A: To construct the Luttinger-Ward functional 



Based on the nonlinear a model approach, v^e analyze an effective Lagrangian Eq. (3). Introducing e'^"" ip with 
the unimodular constraint \ip\-^ = 1, we rewrite Eq. (3) as follows 

C W cl{dr - Hc)Ca + ^K^i + + -^^b{lp*clU + H.C.) 

+ fa{dr - + £/ + A - Wx)/a + ^T^K^j - iai)f„\^ 

Zruf 



2m6 



+^pI - XSN + ^{d^a, - d.a^f, 



(Al) 



where plays the role of an effective chemical potential, imposing the rotor constraint. As a result, three order 
parameters appear to be pb, A, and /^^ beyond the slave-boson mean-field analysis. Introduction of p,^ gives rise to a 
novel energy scale, determining the coherence of the heavy-fermion band. 

One can derive an effective action from our effective field theory Eq. (Al), taking into account quantum corrections 
self-consistently in the Eliashberg approximation, where self-energy corrections are introduced, but vertex corrections 
are neglected. The Eliashberg approximation results in the following effective action 









f d'r' 







3_ 

2m, 



) 5{t - T')S^{r - r') + Se(r - r', r - r')}c4r', r') 



-A/'Se(r - r', T - T')Gc{r' -r,T' - r) 



+ dr J d\J^ dr' J d^r'[f:{r,T)[(^dr-fic + ef + X-^)s{T-T')6^{r-r') + i:f{r-r',T-T')}u{r',T') 

-NJ^fir -r',T- T')Gf{r' ^ r.r' - t)] 

+ [dr j d^rj^r' I d'r'[r{r,T){(^[pb-X/u,]dr-^^d',-p^+PbX-£^^df)5{T-r')5^{r-r') 
+S^(r - r',r - T')}v(r-',T') + Sv,(r - r',T - T')G^{r' -r,T'- r) 

+ jyT j d^rj^T'l d\'[^a^{r,T){(^-^±^)p^J{T-T')5\r-r') + n^4r-r\T-r')}a4r',T') 
+n^,(r -r',T- T')D^,{r' - r, r' - r)] 



-"^ PbJ dr J d\ J dr' J d^r'G^ir - r' ,T - T')Gf{r' - r,T' - T)Gcir - r' ,T - t') 

dr j d^r J\t' J d^r'v^D^^ir - r' ,t - T')vlGf{r' - r,T' - T)Gf{r - r' ,t - t') 
'IJ^^dT J d^rj^^dT'J d^'r'vfD^^ir - r' ,T - T')vtG^{r' - r,T' - T)G^{r - r' ,T - t') 



\2ub 



(A2) 
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where Ec(r — r',T — r'), S/(r — r'.r — r'), S^(r — r',T — r'), and n^,y(r — r',r — r') are self-energy corrections of 
electrons, spinons, ip, and gauge fields, respectively, and Gc{r — r' ,t — r'), Gf{r — r',T — r'), G^(r — r',T — r'), and 
D^uir — r' ,T — t') are their Green's functions. Although vertex corrections are neglected for self-energy calculations, 
such contributions are introduced self-consistently into three coupled equations for order parameters. A way how to 
derive this effective action is shown in Ref. [l^. 

Performing the Fourier transformation and integrating over all field variables, we find the Luttinger-Ward functional 



N f d k c / \ 'I 

jJ2j j^{^ri[-G-\k,tco))+Mk,tco)Geik,tu;)j 

f E / -^{in[~GJ\k,^u;))+^fik,^u;)Gfik,^u;)] 

/ 0^{^^{~G-\q,^n))+J:^iq,^n)G.^iq,^n)} 

^E / ^{ln{-D-\q,^n)) +Uiq,^n)Diq,^n)] 

1 f d^k \ f d^a 
^'P^-pT. J (^^E J J^G.,{q,^n)Gf{kMGc{k + q,^u: + ^n) 



f^E / (0^E / ^Fik,q)Diq,^n)GAkMGf{k + q,^u + ^n) 



4^E / (Sa^E / ^B{q^)D{lMG^[q,^n)G^{q + l^n + iv), (A3) 



where Gc(k,iuj), Gf{k,iuj), G^{k,iui), and D{q,iQ,) are Green's functions of electrons, spinons, phases, and gauge 
fields, respectively, given by 

Gc{k,iuj)^- Gf{k,iuj) 



iuj + lie - - Y.c(k,iuj) iuj + fic - £/ - A - - T,f{k,iuj) 

1 



G^iq,in) - ^ ^ , ^^.^^ 2 



{pb - x/ub){in) - 2^ir - 2^9^ + Pi,- pb^ - if^) 
2— -1- n(9,iii) 



2mb 

D[q,l^) = - o:2^„^ .^. ^ D,,{q,in) = D{q,in)Pf^{q), n,,{q,tn) - n(g, (A4) 



Pij{q) is the projection operator to the transverse component, and 



F{k,q)^-Y, V{{S^J- 
» J = l 
2 



B{q,l) = lj2^t{^ 







f ki+qi/2 

v; = 

ruf 


\ 

P J 


?; - 


%h Qi ~^ 
V — 

rub 



where vf and are velocities of spinons and phases in the i-direction. 



7 



Minimizing the free energy functional with respect to ah seh-energies, we obtain self-consistent Eliashberg equations 



1 f cfik If (Pi 

E^(g,ir!) = / ——Gf[k,iu)Gc{k + q,iuj + in)+pl-Y^ / —---B{q,l)D{lw)G^{q + lin + w), 

iuj ^ ^ iu ^ ^ 

n(g,*r!) = --^ / F{k,q)Gf{k,lLo)Gf{k + q,zuj + zn) + f^-Y, / -^B{q^)G^{q,zn)G4q + 1,1^1 + lu). 

iuj iy 

(A5) 

We note that these equations are essentially the same as those in the quantum critical regime of the Kondo breakdown 
theory [1, |^ , where the self-energy and Green's function of -0 are identified with those of b. 



Appendix B: To evaluate self-energy corrections 

1. Self-energy corrections for the heavy-fermion band 

In order to describe the heavy-fermion band without condensation of tA, we separate fermion self-energy corrections 
as follows 

Y,c{k,ibj) ^ ^c{k,iuj) + l\Y.c{iio), E/(fc, iw) = iw) + AI]/(ia;), (Bl) 

where $c(fc) and iuj) are associated with the formation of the heavy-fermion band, and /S.Tic{iuj) and AI]^(icj) 

are related with non- Fermi liquid physics of such heavy fermions. 

Static contributions of bosons determine the formation of the heavy-fermion band, given by 

= —G^{0,0)Gf{k,iuj) = ———Gf{k,iuj) « — --^-<?/(A:, zw), 

^ ^ \Peff\ ^ lMe//l 

^ ^ \Peff\ ^ Kff\ 



where 



fj.'^f f = fj,^ - PbX~ T,^ {0,0) (B3) 



is an effective chemical potential, essential for coherence. When it touches zero, J should be replaced with KV')!^- 

gc{k,iuj)^- gf{k,iuj)^- p- (B4) 

+ - 2;;r; tuj + fic-€f - X - ^ 

are bare Green's functions. 

Quantum fluctuations of bosons give rise to non-Fermi liquid self-energy corrections of such heavy fermions 



ASc(iw) = — V [ -^-^G^{q,in)Gf{k~q,iu}~in), 

AS/M = -^^E/ i$^G^iq,in)Gaik + q,iu + zn)^^Y. [ i$^F{k,q)Diq,zn)Gf{k,tu) 



(B5) 

where the static component of the ip propagator should not be taken into account. 



8 



For convenience, we also divide the V self-energy as follows 

= S^(0,0) + AS^((7,il7), (B6) 
where the static contribution is introduced into the effective chemical potential, given by 

E^{0,0) = V'pAY1 I 4^Gfi^^'^)Gc{kM^-^^. (B7) 

and the dynamic part results in Landau damping, given by 

1 f (Pk If d^l 

Ai:^{q,in) = V^Pb^Yl / 7^Gf{k,iLo)Gc{k + q,iiO + in) + pl-J2 / ■7^B{q,l)D{l,iu)G^{q + + iu). 

(B8) 

In the next subsection we evaluate this dynamic contribution, where the self-energy correction from gauge fluctuations 
(the second term) will not be taken into account. That contribution is irrelevant because the dynamics is described 
by z = 3, which implies that the effective theory for the xp dynamics lies above the upper critical dimension, resulting 
in the mean-field-like dynamics. See our discussion on this issue in the last section. 



2. To calculate the self-energy 



In order to cover both the quantum critical regime and the incoherent heavy-fermion phase, we introduce the 
following fermion Green's functions 



Gf{i(jj, k) = 



1 



ico — Vp{k — kp) 



Gciico, k) 



TV^Pi./\l^tsj\ 
iu}—Vp.{k—kp) 

1 



AS/(ia;) 



(B9) 



ASy(c)(iw) is the spinon (electron) self-energy due to inelastic scattering with quantum phase fluctuations, where the 

heavy-fermion contribution of <&/((;) (A;, iw) is expressed explicitly. We emphasize that there is T/|/Xgjj|, renormalizing 
/Oft, which reduces the strength of hybridization due to incoherence. We have linearized each bare dispersion, where 
v^p'^^ is the Fermi velocity and k^p'^^ is the Fermi momentum. 

Neglecting non-Fermi liquid parts of self-energy corrections for the time being, we can express Eq. (B9) as follows 



G//c(w,fc) = 



^-vy\k-K,^i'') 



(BIO) 



where Vp^'' and K-'p"^ are renormalized Fermi velocity and renormalized Fermi momentum, respectively, and Z^j^ is 
the wave-function renormalization function. They are given by 



'f/c 



Id 



- (jip + kp"^ 



(k^p-k^) 



c/f 



{vi^)\^.-kf) 



TV^Pt/\ptff\ 

i 7 — 

VpV'}r 



= 1 



TV^Pi,/\ptff\ 



(Bll) 
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Inserting Eq. (B9) with Eq. (BIO) into the ip self-energy, we obtain the following expression 

1 r cPk 

-(.a.vA.,-.^,),o<i±fc|l±ii^) 



(zl7 + V^(/C|,-/C^))2-(aV|,9)2 



+aV^,log(^"^"+">^^>'^>^^ -^""^'^^ )| (B12) 



with a = yl^. Taking a — > 1 with /C^ = /Cj^ for the heavy- fermion band, this expression is simphfied as 

AE,(.17,q) . V^p, 10mn-" + ^^'^+ .f^f + ^1^) + 2V-» ,J-^r , 1(B13) 



If one expands Eq. (B12) in the limit of '^f ^ \^ tj^g typical Landau damping form results 

I]^(zf},g)=7— , (B14) 

q 

originating from particle-hole excitations around the Fermi surface. The damping coefficient is given by 



7 = V'p,^^^^. (B15) 



Then, the ip propagator becomes 



G,iq.^n)^^--^—^—-^, (B16) 



where M^// — f^i> ^ Pb^ — 5]^(0, 0) is an effective chemical potential for phase fluctuations. 

Inserting this boson propagator into self-energy equations for fermions, one can find fermion self-energy corrections. 
Such calculations have been performed in previous studies ^ |^ when the boson dynamics is critical and described 
by z = 3. Since the ip dynamics is also characterized by z = 3 when T > Ti, as discussed in the manuscript, the 
previous results are applied to the present situation directly. Then, we obtain non-Fermi liquid self-energies. 

Appendix C: To derive self-consistent equations for three order parameters of pt, A, and fj,^ 

1. General formulae 

One can find self-consistent equations for order parameters from the Luttinger-Ward free energy functional. An 
essential merit of this approach is that vertex corrections are naturally introduced beyond the Schwinger-Dyson 
equation for an order parameter usually identified with the mean-field equation. 

Minimizing the free energy functional with respect to pf,, we obtain 

where the derivative for each self-energy implies each vertex correction. In particular, we see that ^^^^t^ is the 
vertex correction, performed in the single impurity problem [26| . Such a contribution is expected to modify the 
slave-boson mean-field equation for pb in principle. 
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Minimizing the free energy functional with respect to A, we obtain 

Aro ^V- f dT,c{k,iuj) ^ ,^ . , f <fik / dY,f{k,iu:)\^ ,, . , 



^ E y I- ^i;^" - - dx ) ^""^ + ^ E y (2.)3 = 



where ^^'g^''"'' is the vertex correction beyond the slave-boson mean-field analysis. 
In the same way we find an equation for n^, given by 

N s;-^ f (Pk d'Ec{k,ioj) N ^ f d^k dT,f{k,i 



where ^^°('''^") and ^^/(^''") are identified with vertex corrections. 



(C2) 



2. An equation for yuy, 

We analyze the equation for ix^. Inserting both fermion self-energies associated with the heavy- fermion band and 
boson self-energy into Eq. (C3), we obtain 

f d^k r d fV^pt T . f d^k r d fV^pt T . . 



We rearrange this equation as follows 

I f d^q -o^ ^ 1 ^ /■ '^''^ / ^ (V^PblCf\ d / If^K^ 

1 + ^Ey + ^Ey (^la^i^li;^ v^j - ajr.l^YJr^^'''''^ 

(C5) 

This expression is quite interesting in the respect that the first two terms in the left-hand-side correspond to the 

Schwinger-Dyson equation resulting from 1 while other contributions originate from vertex corrections. 

Keeping /x^-derivative terms only when they depend on /x^ explicitly as the lowest-order approximation, we reach 
the following expression 



1-^^{;^E J ^9f{k,icv)Gc{k,iuj) + ^Yl J (||35-(^'*'^)<^/(^'^^)} = -^E J ^^^^"l'^^^- 

(C6) 



It is clear that fermion contributions are related with vertex corrections. We will see that this correction plays an 
important role for self-consistency, which cancels other quantum corrections. 
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3. An equation for pb 

Inserting both fermion self-energies associated with the heavy- fermion band and boson self-energy into Eq. (CI), 
we obtain 

f d^k c d fV^pb T . ^^ . TV ^ f d^k c d (V^pb T . x-, 

+5i:/(4{'"-4-* + 5^(l^Tf)-5;;;(-'V)K<'''") = °- (• 



Performing derivatives for p^, we reach the following expression 

q^Y^ ( 1 ^ [ d^k 1 . f d^k "1 

^bPb-—--^{-^Yl J j2^9f{k,iuj)Gcik,iw) + - Y, J (2^5c(A;,ia;)G/(fc,ia;)| 



+ 



Pb /^ V-^ K.p \ ( I sr^ f d'^k . . ^ ^ sr^ f di^k ,, . . o 



in ^ ^ in ^ ' ^ in 

Inserting Eq. (C6) into the above equation, we obtain 

'1 — d^k . . . I ^ f d^k 



j^Y ( 1 , f d k 1 ^ f d k "1 

y j^g}{k,ioj)Gc{kM + ^Y. J ■^9cik,iw)Gf{k,iuj)j 

f'^eff ^ iu) ' iuj ^ ^ 



+A 



^^{;^E / -^9f{kMGc{k,iuj) + \Y. I 4i39cik,iu;)Gf{k,iuj)}] 



1 



Surprisingly, quantum corrections in the tp sector cancels those in the fermion part, simplifying the above expression 
as follows 



TV (1^ f d-'k n ■ n ■ ^ f d:'k . . ,\ , V K."/ 



This cancellation confirms the validity of our approximation in Eq. (C6). 

Neglecting the right-hand-side because we approximate the gauge propagator as Eq. (B16), where the linear time- 
derivative is not introduced, we reach the following expression 



A-^^ = ^{^E y J^9AkMGe{k,^.) + -Y: J j^9c{kMGf{kM}, (Cll) 



essentially the same structure as that of the slave-boson mean-field theory except for the second term in the left-hand- 
side. 
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4. An equation for A 



Inserting both fermion self-energies associated with the heavy- fermion band and boson self-energy into Eq. (C2), 
we obtain 

N ^ r cPk c d fV^pb T . ^^^^ „ . ^ 

v^y (2^{i-aA(^^^^(^'^^))}^^(^'^^) 

Performing A-derivatives for terms depending on A explicitly as the lowest-order approximation, we obtain 

J J^9f{k,iuj)Gc{k,iu;) + -J2 J j^gc{k,iu;)Gfik,zuj)\ 

f d^k ^ . , r TV^pbfl^ f d^k n ■ r, ■ , 1 V- /■ . „ . 

+ / 7^'^/('''''^^ + ^n^~Tv^i;s2^ / 7^9f{k,iuj)Gc{k,iu) + / — -3 5c(A;,^a;)G/(fc,^a;)| 

1./,) ^ ' r*f> f f ifii \ / '"^ A,,, v V / 



Also, quantum corrections in the ip sector cancels those in the fermion part, recovering the constraint equation in the 
slave-boson mean-field analysis 

N I' fP k 

when the first term in the right-hand-side of Eq. (C13) is neglected. This treatment is consistent with the 'ip Green's 
function, where the X/ub term in the linear time-derivative is not considered. 



Appendix D: To solve self-consistent equations for order parameters 



Three coupled self-consistent equations are given by 



/ ^9fi''^i'')<^c{k,iuj) + J -^9c{k,iuj)Gf{k,iuj)} = ~Y1 J (^'^'/'(a. 
^-^^ = ^{^E/(|^5/(A'-)Gc(A'-) + ^E/ ^9cik,^u.)GAk,^u:)}, 



p^ + ^E 



il^k 

j^Gf{kM = NS 



(Dl) 



beyond the mean-field analysis, where quantiim corrections are introduced self-consistently. 

Inserting both renormalized lieavy-fcrmion Green's functions and renormalized ^ propagator into Eq. (Dl), we 
reach the final formulae for two energy scales 



1- 



Pb 



d^q 



I// I 



-I- 



^e// 



\ £_ _ o 

27r2 V% l-.V- 



e//l 



dk ^^f{E+{k))- f{E^{k)) 



27r2 



E+{k)-E_{k) 



Pb + 2 



E+{k)-E_{k) 



E+{k)-E_{k) 



(D2) 
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where the effective chemical potential and the renormalized heavy-fermion band are 



27r2 VS 



E±{k) = l{ef{k)+s,{k)) ± y{ef{k)-e,{k)y+4TV^pb/\ntffl 



respectively. The damping coefficient is given by 



47rV^V^ 



where renormahzed Fermi momentum and renormahzed velocity are 



(D3) 



(D4) 



— K-p — - (jcp + kp^ "'"2^ (^^ ^"^^ 



VttV 



pvp 



'Sic 



(vfy{Kg^-kfy\ 



z-j = 1 



TV'pb/\^^tff\ 



(v^j^y{ici^^-kry 



(D5) 



/(e) is the Fermi-Dirac distribution function. Then, all quantities are defined, where nic, mf, e/, and /ic ai'e only 
parameters to define our system. 

When vertex corrections are neglected in these equations, we obtain 



1 



1 



ry2 



00 

dk ,,f{E+{k))-f{E^{k)) 



P» + '^l^k^ 



E+{k)-E_{k) ' 



E+{k) - E-{k) 



E+{k)-E-{k) 



(D6) 



where the second and third equations recover the slave-boson mean-field equations, replacing T/\p'^jj\ with KV")!^ in 
the second equation. 

One may consider that the first equation of the rotor constraint incorporates quantum corrections fully self- 
consistently because the tp self-energy correction is introduced. On the other hand, the fermion self-energy for 
non- Fermi liquid physics is not taken into account. In our opinion introduction of such self-energy corrections will 
not change the present picture, in spite of modifying it only quantitatively. 



Appendix E: Numerical analysis 
The higher energy scale T* is determined from Pb{T*) = 0. Then, self-consistent equations are reduced to 



i = T*T L 



d-^l }^ * _ T/2 kpkp 

in V 7 g p-gff VpVp 



T*V^ f dk f*(efik))-ne,ik)) 



V'kf_ T*V' fd^ 

27r2^;^ lM^;/l/27r2 ef{k)-e,ik) ' 

00 

/r]k 
^k'n<^Ak)), (El) 
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where the conduction band is decoupled from the spinon band. Notice that the boson band becomes flat, resulting in 
incoherence as long as —fJ-tjf > 0- 

Solving the third equation, we obtain A as a function of T*. Inserting the A into the second equation, we find 
^*/lMe//l ^ function of both T* and V. Inserting this function into the first equation, we obtain an equation, 
representing the relation between T* and V, where the following cutoff scheme is used. 



dq 



1 



J^n JO 



1 



(E2) 



As a result, we find the T*{V) line in the phase diagram of Fig. 1. 

It is subtle to determine the Fermi-liquid temperature. It is identified with the condensation temperature of thus 
given by f^^j^fiTpL) = 0. Since Tpi^/\ii^jj:{Tpl)\ diverges, self-consistent equations are not well defined. It is natural 

to replace T/|/i^^j| with KV-")!^- Taking fjifj^ = with {tp) = 0, we reach the following equations to determine TpL, 



E 

in 



7 = V^pb 



\FL 



27r2 VS, 



2V^ 



dk 
2^ 



E+{k) - E_{k) 



Pb 



dk 
2^ 



.{k-k'p)^E^{k) 



E+{k) - v^p{k - k^p) 



E+{k) - E_{k) 



E+{k)-E_{k) 



(E3) 



where 



E±{k) = hef{k)+e,{k)) ± 




iV'pt)vf 



)( 



+ 4: 



V^Pb 



F^F 



= 1 



V'pb 



K,p — Kp 



2 • 



(E4) 



One may criticize that Eq. (E4) are nothing but those of the heavy-fermion state. Since T/|/i^^^| is replaced with 
one can claim that V^pt should be modified into V'^pb\{'ip)\'^ ■ However, this substitution gives rise to a serious 
problem. When (■0) vanishes at T = TpL, the resulting state does not have the heavy-fermion band. This conclusion 
is in contrast with the existence of T*, where the heavy-fermion band results, but its coherence is not achieved yet. 
In this respect we perform the above approximation, where we resort to the "bare" heavy-fermion band. Although to 
use this band structure overestimates TpL, we believe that the final conclusion will not change at least qualitatively. 

It is straightforward to solve Eq. (E3) because the last two equations are decoupled from the first equation. The 
last two equations are nothing but the slave-boson mean-field equations. Solving these coupled equations, we obtain 
Xfl and pb as a function of both TpL and V. Inserting these functions into the first equation, we can determine 
TfUV) in Fig. 1. 



Appendix F: Remarks 



1. On decomposition for the RKKY spin-exchange term 



One may criticize the way how to take the RKKY spin-exchange interaction in terms of spinous from Eq. (1) to 
Eq. (2). A systematic description is to introduce the Sp(2N) representation, which allows us to construct the spin op- 
erator in the large-N case [Ying Ran and Xiao-gang Wen, arXiv:cond-mat/0609620 . Decomposing the spin-exchange 
interaction into both spin-singlet and spin-triplet channels and performing the Hubbard-Stratonovich transformation 
for both exchange hopping and pairing interactions in each spin channel, one can find an effective Hamiltonian for 
paramagnetic Mott insulating phases (or spin liquids) [Ryuichi Shindou and Tsutomu Momoi, Phys. Rev. B 80, 
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064410 (2009)]. All these possible mean-field phases for spin liquids can be classified in terms of the so called projec- 
tive symmetry group [Xiao-Gang Wen, Phys. Rev. B 65, 165113 (2002)]. Performing the classification to distinguish 
various spin liquid states, one takes into account the Gutzwiller projection for such mean-field phases or integrates 
over gauge fluctuations accurately, and finds actual paramagnetic ground states in each parameter regime. This is 
the full procedure along this description. 

However, this full procedure gives rise to much complication in the analysis. Instead, we usually resort to experi- 
ments, helping us neglect several irrelevant interactions in the system, which will be determined by the Fermi-surface 
structure basically. The Kondo breakdown scenario [I. Paul, C. Pepin, and M. R. Norman, Phys. Rev. Lett. 98, 
026402 (2007); C. Pepin, Phys. Rev. Lett. 98, 206401 (2007)] has been applied to quantum critical physics in 
YbRh2Si2, and successfully explained thermodynamics [K.-S. Kim, A. Benlagra, and C. Pepin, Phys. Rev. Lett. 
101, 246403 (2008)], both electrical and thermal transport [K.-S. Kim and C. Pepin, Phys. Rev. Lett. 102, 156404 
(2009)], and the Seebeck coefficient [Ki-Seok Kim and C. Pepin, Phys. Rev. B 81, 205108 (2010); K.-S. Kim and 
C. Pepin, Phys. Rev. B 83, 073104 (2011)] although this theory underestimates antiferromagnetic correlations. Re- 
call that the Kondo breakdown theory is based on the uniform spin-liquid ansatz for the antiferromagnetic phase, 
where quantum critical physics is described by Fermi surface (Kondo) fiuctuations, expected to be applicable to the 
finite-temperature regime at least. This quantum critical physics is described by the Eliashberg approximation, where 
quantum corrections from conduction electrons, spinous, holons, and even gauge fluctuations are incorporated self- 
consistently. Maybe, the renormalization group analysis will make our understanding of the present subject deepen. 
But, we try to justify our effective Hamiltonian within the phenomenological background. 



2. On phase fluctuations in the single impurity problem 



It is necessary to review the 1/N correction for the Kondo problem in more detail. In particular, one may criticize 
our description on the linearization for the phase factor, i.e., e'^'' ~ 1 + iOi, because this procedure seems to break the 
gauge invariance, which may cause unphysical divergences. In the "Higgs" phase such phase fluctuations should be 
eaten by gauge fluctuations, giving rise to the fact that all field variables in the resulting effective Hamiltonian are 
given by gauge singlets, where the remaining gauge fluctuations (in the unitary gauge) are gapped, and do not affect 
the low energy physics. But, this should be regarded as just one way to describe such phase fluctuations. 

Actually, the linearization was proposed in N. Read, J. Phys. C: Solid State Phys. 18, 2651 (1985). Although this 
seminal paper focuses on fluctuation corrections (1/N) in the complex boson representation instead of the angular 
or polar coordinate representation, the angular representation is also discussed (section 6). As well discussed in the 
manuscript, such quantum corrections give rise to the Zog-divergence for the expectation value of holon condensation 
in the complex boson representation, which implies that the condensation does not appear, where the expectation 
value vanishes as the cutoff goes to inflnity. Actually, the holon propagator exhibits the power-law dependence in 
time at long time scales instead of a constant. 

Before we turn to the angular representation, we would like to emphasize that an important point is how quantum 
corrections are incorporated self-consistently well, keeping the Ward identity or conservation law in any representa- 
tions. Even if the perturbation is used in the angular representation, the final result will be physically valid as long 
as the Ward identity is satisfied. The linearization gives rise to new type of a vertex, resulting in the additional 
log-divergent contribution to the spinon self-energy. This ^0(7-divergence was argued to cancel the /og-divergence 
in the expectation value of the holon condensation, where the holon-condensation amplitude becomes finite. This 
means that the holon condensation (complex number) vanishes due to transverse (phase) fiuctuations while the holon- 
condensation amplitude (real number) remains finite. In this respect our present work can be regarded to generalize 
the single-impurity problem to the impurity-lattice problem. Indeed, our methodology is to generalize that of the 
Kondo problem [N. Read, J. Phys. C: Solid State Phys. 18, 2651 (1985)] to the lattice problem. 

The reason why we took this way is that we want to describe the crossover regime out of the Higgs phase, 
where strong phase fluctuations make the unitary gauge useless because singular configurations for phase fluctu- 
ations, for example, vortices, should be taken into account in that description [F. S. Nogueira and H. Kleinert, 
'arXiv:cond-mat / 0303485, the World Scientific review volume "Order, Disorder, and Criticality" , Edited by Y. Holo- 
vatch]. In this case the unitary gauge is not an easy way to describe the crossover regime. Instead, we resort to 
the non-linear cr-model description, allowing us to introduce higher order quantum corrections and to describe the 
crossover regime. 
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3. On the holon self-interaction term 

It is necessary to discuss the physical origin of the holon self-interaction term, phenomenologically introduced in the 

study. Before this discussion, wc would like to mention that this term does not play any important role in low energy 
critical dynamics because the Landau damping term is mostly relevant, which results from dynamics of fermions near 
the Fermi surface. 

Microscopically, the self-interaction term can originate from both the J-tcrm, i.e., J'^^jiSi ■ Sj — jUirij), where 
the density-interaction term with the density operator rii of localized electrons is not shown explicitly in the present 
manuscript, and the dynamically generated self- interaction term resulting from quantum fluctuations of fermions, 
where integrating over fermions and expanding the resulting logarithm up to the fourth order for holon excitations 
give rise to such a term. The second effect is discussed in I. Paul, C. Pepin, M. R. Norman, Phys. Rev. B 78, 035109 
(2008); C. Pepin, Phys. Rev. B 77, 245129 (2008). Since this effect appears from the short-distance scale, such 
physics is not universal, irrelevant at low energies. 

4. On fluctuations 

One may criticize our approximation to neglect /i^, fluctuations for the ijj dynamics. Actually, such fluctuations 
play an important role in the nonlinear a model description of the XY-type model, where the boson dynamics is 
described by 2; = 1, i.e., the relativistic dispersion [T. Senthil, Phys. Rev. B 78, 045109 (2008)]. However, our phase- 
fluctuation dynamics is governed by z = 3 due to Fermi-surface fluctuations, as emphasized before. As a result, the 
boson dynamics can be understood within the "mean-field" approximation, which neglects /i^ fluctuations, because 
the effective field theory for the boson dynamics lives above the upper critical dimension, i.e., d + z > dc = 4:. In the 
same way gauge fluctuations do not affect or alter the dynamics of boson excitations although they axe responsible 
for anomalous self-energy corrections to spinon excitations. Although the situation is not completely the same as 
our case, wc would like to refer to T. Senthil, Phys. Rev. B 78, 045109 (2008) for irrelevance of gauge-fluctuation 
corrections to the holon (phase) self-energy. 

Even if one focuses on the case {z = 1) when fluctuations are relevant, such fluctuations cannot erase the 
existence of the incoherent heavy-fermion phase. They can modify the critical physics in the preformed heavy-fermion 
phase, where some critical exponents may be changed. 



